Peak power reduction device in communication system

ABSTRACT

A peak power reduction device includes a unit for dividing digital information to be transmitted into a plurality of streams; a unit for selecting the modulation level of the streams and distribution of transmission power according to a transmission state; a unit for performing singular-value decomposition on the transmission path characteristic of a streams and precoding the resultant data by a right singular value matrix; a unit for performing complex mapping on the subcarrier of a stream according to the modulation level; a unit for converting a complex mapping signal into a time domain signal; a first unit for storing a conversion result as a time domain signal; a second unit for calculating a peak time signal exceeding a predetermined threshold value from a peak value, if any exists, that exceeds a prescribed threshold value for the amplitude of the time domain signal.

FIELD OF THE INVENTION

The present invention relates to a method of reducing peak power versus average power in eigenmode MIMO transmission for transmitting a plurality of streams of a communication system and, more particularly, to a device for peak power reduction in multicarrier transmission such as an OFDM scheme or the like.

BACKGROUND OF THE INVENTION

A MIMO (Multiple Input Multiple Output) technique that spatially multiplexes and transmits a plurality of streams using a plurality of transmitting/receiving antennas has been put to practical use. Transmission capacity can be increased by this MIMO transmission. As one of MIMO transmission schemes, there is available an eigenmode MIMO transmission scheme that can use, to the utmost, the potential transmission capacity of a transmission path when the state of a transmission path can be grasped on the transmission side.

In the eigenmode MIMO transmission scheme, assuming the transmission path characteristic as H as shown in Expression 1, the transmission path characteristic is singular-value-decomposed as preprocessing of eigenmode transmission.

H=UΣV ^(H)  (Expression 1)

(H on the upper right side of V indicates complex conjugate transpose)

In Expression 1, U denotes a left singular value matrix, V denotes a right singular value matrix, and Σ is a matrix indicating an eigenvalue (√λ_(k)) whose diagonal component is H as shown in Expression 2. In addition, the elements of each matrix are complex numbers. In the following description, the elements of a matrix and a vector are assumed to be complex numbers unless otherwise stated.

Σ=diag{√{square root over (λ₁)}√{square root over (λ₂)} . . . √{square root over (λ_(L))}}  (Expression 2)

where λ₁>λ₂> . . . >λ_(L).

In eigenmode MIMO transmission, L stream vectors S=[S1 S2 . . . SL]^(T) are used as streams for performing information transmission. Information is transmitted by placing the information on the amplitude and phase of the stream vectors S. After multiplying the stream vectors S by the right singular value matrix V, a transmission signal vector X is generated. Furthermore, the transmission power of each stream is appropriately distributed by using the principle called a water injection theorem, where by the increase in the transmission capacity can be expected.

Assuming that the power distribution matrix is P=diag{√{square root over (p₁)} √{square root over (p₂)} . . . √{square root over (p_(n-1))}}, the transmission signal vector X is represented by Expression 3.

X=VPS  (Expression 3)

The processing shown in Expression 3 is referred to as precoding processing on the stream vectors S. The transmission signal vector X after precoding is transmitted from each transmission antenna.

The reception signal vector Y is a signal obtained by adding a thermal noise vector N (where a variance is σ2) after the transmission signal vector X passes through the transmission path matrix H.

Y=HX+N  (Expression 4)

Expression 1 and Expression 3 are substituted into Expression 4 to obtain Expression 5.

Y=(UΣV ^(H))(VPS)+N  (Expression 5)

The transmission stream vector S indicated in Expression 6 is estimated by subjecting the reception signal vector Y represented by Expression 5 to spatial filter processing for decomposing the reception signal vector Y for each stream. The spatial filter processing can be realized by multiplying the reception signal vector Y by a complex conjugate transpose matrix UH of the left singular value matrix.

Ŝ=U ^(H) Y

Ŝ=U ^(H)(UΣPS+N)

Ŝ=ΣPS+U ^(H) N  (Expression 6)

In Expression 6, the first term ΣPS on the right side denotes a signal component, and the second term U^(H)N denotes a noise component. The details of the signal component of the first term are shown in Expression 7.

$\begin{matrix} {{\Sigma \; {PS}} = {\begin{bmatrix} \sqrt{\lambda_{1}} & 0 & 0 & 0 \\ 0 & \sqrt{\lambda_{2}} & 0 & 0 \\ \vdots & \vdots & \vdots & \vdots \\ \vdots & \vdots & \vdots & \vdots \\ 0 & 0 & 0 & \sqrt{\lambda_{L}} \end{bmatrix}{\quad{\left\lbrack \begin{matrix} \sqrt{P_{1}} & 0 & 0 & 0 \\ 0 & \sqrt{P_{2}} & 0 & 0 \\ \vdots & \vdots & \vdots & \vdots \\ \vdots & \vdots & \vdots & \vdots \\ 0 & 0 & 0 & \sqrt{P_{L}} \end{matrix} \right\rbrack {\quad{{\begin{bmatrix} s_{1} \\ s_{2} \\ \vdots \\ s_{L} \end{bmatrix}\mspace{20mu} \Sigma \; {PS}} = \begin{bmatrix} {\sqrt{\lambda_{1}P_{1}}s_{1}} \\ {\sqrt{\lambda_{2}P_{2}}s_{2}} \\ \vdots \\ {\sqrt{\lambda_{n - 1}P_{n - 1}}s_{n - 1}} \end{bmatrix}}}}}}} & \left( {{Expression}\mspace{14mu} 7} \right) \end{matrix}$

Similarly, the details of the noise component of the second term are shown in Expression 8.

$\begin{matrix} {{U^{H}N} = {U^{H}\begin{bmatrix} n_{1} \\ n_{2} \\ \vdots \\ n_{n - 1} \end{bmatrix}}} & \left( {{Expression}\mspace{14mu} 8} \right) \end{matrix}$

In Expression 8, [n1 n2 . . . nn−]⁷ is a thermal noise exhibiting a Gaussian distribution, the average value of each element is 0, and the variance is σ2. Since UH is a unitary matrix, there is no such thing as emphasizing this noise level.

The SN ratio of each stream is calculated using Expression 7 and Expression 8. Assuming that the SN ratio of the k^(th) stream is SNRk, Expression 9 is obtained.

$\begin{matrix} {{{SNR}_{k} = \frac{\lambda_{k}P_{k}{E\left\lbrack {s_{k}}^{2} \right\rbrack}}{\sigma^{2}}}{{SNR}_{k} = {\lambda_{k}P_{k}\frac{P_{s}}{\sigma^{2}}}}{{{where}\mspace{14mu} P_{s}} = {{E\left\lbrack {s_{k}}^{2} \right\rbrack}.}}} & \left( {{Expression}\mspace{14mu} 9} \right) \end{matrix}$

As shown in Expression 9, if appropriate precoding processing is performed, the SN ratio of each stream becomes a value corresponding to the eigenvalue of the transmission path and the power distribution matrix. Although not described in detail, according to the theory of the water injection theorem of the eigenmode MIMO transmission scheme, it is possible to transmit a large amount of information by performing large power distribution √P_(k) for a stream having a large eigenvalue √λ_(k).

Therefore, by allocating a modulation scheme with a large multi-value number (for example, 256 QAM (Quadrature Amplitude Modulation), etc.), it is possible to make maximum use of the potential transmission capacity of the transmission path. Conversely, since it is impossible to transmit a lot of information with a stream having a small eigenvalue, a modulation scheme with a small multi-value number (for example, QPSK (Quadrature Phase Shift Keying), etc.) is allocated.

By the eigenmode transmission MIMO scheme described above, it is possible to perform transmission with the maximum transmission capacity for the transmission path.

In addition, when the eigenmode MIMO transmission is applied to a wide frequency bandwidth, the characteristics of the transmission path often become different for each frequency ω due to the influence of multipath or the like. Therefore, the transmission path characteristic matrix H(ω) is represented by a function of the frequency ω. In an environment having such frequency characteristics, it is possible to optimize the transmission capacity by appropriately performing the above-described stream distribution for each frequency (ω). In order to perform optimization on a frequency-by-frequency basis, an OFDM (Orthogonal Frequency Division Multiplexing) modulation scheme that performs transmission while making narrowband multiple subcarriers orthogonal to each other is used, whereby stream distribution suitable for the transmission path H (ω can be performed for each subcarrier (or on a unit of a plurality of subcarriers).

By the way, in the OFDM scheme, when the phases of a plurality of subcarriers coincide with each other, there is a tendency that a large amplitude peak (instantaneous peak) appears.

On the other hand, in the case of a general power amplifier, when a signal with a large amplitude is inputted and output power comes closes to a saturated state, nonlinear distortion occurs. This also applies to a power amplifier used for amplifying a carrier on the transmission side of a transmission system. As in the case of the above-described OFDM scheme, if a large amplitude peak is included in the time waveform of an input signal, nonlinear distortion of a power amplifier occurs. Due to this distortion, the EVM (Error Vector Magnitude: modulation accuracy) of all streams is deteriorated, and the bit error rate is increased. Furthermore, due to this nonlinear distortion, there is also posed a problem that power leakage out of a band occurs.

Although the problem attributable to the nonlinear distortion can be solved by increasing the back-off of the power amplifier, there is a tradeoff relationship that the efficiency of the power amplifier decreases. Thus, as a technique for reducing the instantaneous peak, a technique has been proposed in which filtering is performed on a peak exceeding a predetermined level of threshold value, thereby reducing the peak while suppressing out-of-band leakage (see, e.g., Patent Document 1, etc.).

PRIOR ART DOCUMENT

Patent Document 1: Published Japanese Translation of PCT International Application No. 2009-543434

SUMMARY OF THE INVENTION

In the above-described conventional technique, it is necessary to insert a signal for peak reduction into the entire radio signal band. Therefore, EVM deterioration is inevitable. In addition, as described above, when eigenmode MIMO transmission is performed, a high multi-value number is allocated to a stream which can obtain a large eigenvalue. Thus, a bit error rate increases even with a slight EVM deterioration. On the other hand, a low multi-value number is allocated to a stream of small eigenvalue. Therefore, even if EVM deterioration occurs, a bit error rate does not increase greatly.

It is therefore an object of the present invention to control the EVM deterioration for each stream in accordance with the modulation multi-value number, suppress the EVM deterioration of a stream with a large eigenvalue, and conversely allow the EVM deterioration of a stream with a small eigenvalue, so that the peak power reduction can be achieved without deteriorating the bit error rate of the entire eigenmode MIMO transmission.

In accordance with an aspect, there is provided a peak power reduction device including: a unit for dividing digital information to be transmitted into a plurality of streams; a unit for selecting a modulation multi-value number of the plurality of streams and distribution of transmission power according to a transmission state; a unit for singular-value-decomposing transmission path characteristics of the streams and precoding the transmission path characteristics of the streams by a right singular value matrix; a unit for complex-mapping subcarriers of the streams according to the modulation multi-value number; a unit for converting a complex mapping signal into a time domain signal; a first unit for storing a conversion result as the time domain signal; a second unit for, when a peak value exceeding a predetermined threshold exists with respect to an amplitude of the time domain signal, calculating a peak time signal exceeding a predetermined threshold value from the peak value; a third unit for converting the peak time signal into a frequency domain signal and calculating a converted frequency domain peak signal; a fourth unit for calculating a peak stream signal by multiplying the frequency domain peak signal by a complex conjugate transpose matrix of a right singular value matrix; and a fifth unit for limiting an SN ratio of a peak stream signal level corresponding to the modulation multi-value number of the streams with respect to the peak stream signal obtained by the fourth unit to a predetermined deterioration amount or less.

The peak power reduction device may further include: a sixth unit for calculating a peak time signal whose level is restricted by converting the level-restricted peak stream signal to a time domain; and a seventh unit for reducing a peak component by subtracting the time domain peak signal obtained by the sixth unit from the time domain signal obtained by the first unit.

In the peak power reduction device, a modulation error ratio for each stream may be controlled according to the modulation multi-value number allocated to the streams by repeatedly operating the second unit to the seventh unit to calculate a peak-reduced time domain signal multiple times.

Effect of the Invention

According to the present invention, it is possible to reduce the peak power without deteriorating the bit error rate of the entire eigenmode MIMO transmission.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a first example of a peak power reduction device according to an embodiment of the present invention.

FIG. 2 is a diagram for explaining a configuration of stream level control integration of a peak power reduction device according to an embodiment of the present invention.

FIG. 3 is a diagram for explaining a subcarrier configuration.

FIG. 4 is a conceptual diagram of a peak signal.

FIG. 5 is a diagram for explaining conversion between a stream region signal and a transmission frequency domain signal.

FIG. 6 is a diagram for explaining a peak stream vector and a band-limited peak stream vector.

FIG. 7 is a diagram for explaining the equivalence between calculation of a stream region and calculation of a time domain.

FIG. 8 is a diagram for explaining Expression 28 and an example of an allowable EVM deterioration amount.

DETAILED DESCRIPTION OF THE EMBODIMENTS

Hereinafter, a peak power reduction device according to the present invention will be described in detail with reference to the embodiments shown in the drawings.

FIG. 1 shows a first example of a peak power reduction device according to an embodiment of the present invention. The peak power reduction device is provided on the input side of a power amplifier in order to reduce the peak appearing in an input signal of the carrier amplifying power amplifier existing on the transmission side of, for example, an OFDM wireless transmission system.

In addition, as shown in FIG. 3, the number of points of FFT (Fast Fourier Transform) used for OFDM transmission is assumed to be K_(FFT), and the number of subcarriers used is assumed to be K. At this time, it is assumed that the subcarriers in a transmission band are composed of a DATA carrier (information carrier) for transmitting information and a NULL carrier (out-of-band carrier). The DATA carrier is mapped by a modulation scheme such as QAM (Quadrature Amplitude Modulation) or the like and is used for information transmission.

The present invention will be described with reference to an OFDM signal for a radio transmission system that performs eigenmode MIMO transmission using the above-described OFDM transmission scheme. However, it may be possible to use a multicarrier transmission scheme similar to OFDM.

In the eigenmode MIMO transmission scheme, the precoding matrix shown in Expression 3 is required on the transmission side. Before explaining the peak reduction which is the main object of the present invention, reference is made to a method of obtaining the precoding matrix. The method of obtaining the precoding matrix includes two major methods.

The first method is a method of performing a singular value decomposition calculation on the transmission path characteristic matrix H(ω) obtained on the reception side and sending the precoding matrix and the eigenvalue matrix back to the transmission side. This is a method of sending the precoding matrix close to its characteristics back to the transmission side as a code book with respect to the precoding matrix obtained on the reception side.

The second method is a method of acquiring the precoding matrix by singular-value-decomposing the transmission path characteristic matrix H(ω) with the singular value decomposition 101 on the transmission side as shown in Expression 2. In the second method, it is indispensable that the transmission path characteristic matrix H(ω) is obtained on the transmission side. A plurality of realization methods is proposed as the method of obtaining the transmission path characteristic matrix H(ω) on the transmission side. For example, there are available a method of directly transmitting the transmission path characteristic matrix H(ω) obtained on the reception side to the transmission side and a method of obtaining the transmission path characteristic matrix H(ω) using the duality of the transmission path when the time division multiplexing scheme is used. In the present invention, the method of obtaining the transmission path characteristic matrix H(ω) will not be particularly mentioned. The following description will be made on the premise that the transmission path characteristic matrix H(ω) has been obtained.

The transmission path characteristic matrix H(ω) is a matrix of NT rows and NR columns, where NT and NR indicate the number of transmitting antennas and the number of receiving antennas. In this regard, the element of (i, j) of H(ω) is the transmission path characteristic between the i^(th) transmitting antenna and the j^(th) receiving antenna.

The transmission path characteristic matrix H(ω) is inputted to the singular value decomposition 101. In singular value decomposition 101, the transmission path characteristic matrix H(ω) is decomposed into a left singular value matrix (U(ω)), a right singular value matrix (V(ω)) and an eigenvalue matrix (Σ(ω)) of H(ω) as shown in Expression 1. In this regard, the eigenvalue matrix Σ(ω) is represented by Expression 10.

Σ(ω)=diag{√{square root over (λ₁(ω))}√{square root over (λ₂(ω))} . . . √{square root over (λ_(L)(ω))}}  (Expression 10)

The eigenvalue matrix (Σ(ω)) of the transmission path characteristic matrix H(ω) is inputted to the stream multi-value number control (102). On the other hand, the digital information C(ω) to be transmitted is inputted to the mapping (103). The digital information C(ω) is transmitted by placing the digital information C(ω) on the amplitude or phase of a constellation such as QAM (Quadrature Amplitude Modulation) or PSK (Phase Shift Keying). This process is called mapping. As the multi-value number for mapping increases, the information that can be transmitted becomes larger. However, the distance between signal points becomes narrower. Therefore, errors due to a noise are likely to occur. In performing this mapping process, the multi-value number is selected based on the multi-value number selection information from the stream multi-value number control (102). As described above, in the eigenmode MIMO transmission, as the eigenvalue (λ_(i)) of the transmission path becomes larger, the power of the transmission signal is also allocated larger. The transmission capacity can be increased by allocating a large scale multi-value number. Accordingly, the stream multi-value number control (102) allocates a multi-value number of each stream depending on the size of each element √λ₁(ω) √λ₂(ω) . . . √λ_(L)(ω) of the eigenvalue matrix (Σ(ω)) from the singular value decomposition (101). The stream multi-value number control (102) performs the power distribution for each stream at the same time and outputs the power distribution matrix P(ω) to the power distribution matrix (106).

As an example of the multi-value number allocation, 64 QAM (6 bit transmission) is allocated to the first stream having the largest eigenvalue, 16 QAM (4 bit transmission) is allocated to the second stream having the second largest eigenvalue, and QPSK (2 bit transmission) is allocated to the third stream having the smallest eigenvalue.

The mapped stream vector S(ω)=[s₁(ω) s₂(ω) . . . s_(L)(ω)]^(T), to which a proper multi-value number is allocated, is inputted to the 0 padding 104.

The 0 padding 104 is a process provided for performing oversampling in the inverse Fourier transform 108. The oversampling is performed for the number of FFT points, which is calculated by the inverse Fourier transform 108, at a rate of an integer multiple of (KFFT) (generally, a power of 2). For example, in the inverse Fourier transform with KFFT=1024 points, the number of oversampling FFT points is calculated with 4 KFFT=4096 points which is four times. Although the necessity of 0 padding process will be described later, it is desirable that the oversample is four or more times of the number of FFT points. At this time, in the 0 padding 104, 0 is additionally added just as much as (4−1) KFFT=3 KFFT.

The output signal of the 0 padding 104 is inputted to the selection 105. In the process of the present invention, in order to reduce the peak of the OFDM signal, it is necessary to repeat, plural times, a series of processes from the power distribution matrix 106 to the stream level control integration 118 shown in FIG. 1. First, the first repetition process will be described and, then, the second and subsequent processes will be described.

In the selection 105, the signal from the 0 padding 104 and the signal from the stream level control integration 118 are selected and outputted to the power distribution matrix 106. In the selection 105, the signal from the 0 padding 104 is selected and outputted in the first repetition process. In the second and subsequent processes, the signal from the stream level control integration 118 is selected and outputted.

In the power distribution matrix 106, the stream vector S(ω) is multiplied by the power distribution matrix P(ω). The power distribution matrix P(ω) to be multiplied at this time is based on the result calculated by the stream multi-value number control 102.

Thereafter, in the right singular value matrix 107, the right singular matrix V(ω) from the singular value decomposition 101 is multiplied to obtain a transmission signal vector X(ω).

These processes are similar to Expression 3 and are re-listed in Expression 11 as a function of the frequency ω.

X(ω)=V(ω)P(ω)S(ω)  (Expression 11)

The transmission signal vector X(ω) is inputted to the inverse Fourier transform 108 and is transformed from the signal in the frequency domain (ω) into the signal in the time domain (t). At this time, as described in the 0 padding 104, while conducting oversampling, the inverse Fourier transformation is performed with a point number of an integer multiple of KFFT. In order to implement the reduction processing of the peak component of the OFDM signal with high accuracy, it is desirable to perform the reduction processing in the region of an analog signal. However, since the signal processing in the analog domain is difficult, the processing is performed with the digital signal after oversampling in order to make the analog signal as close as possible.

This inverse Fourier transform processing is shown in Expression 12.

x ⁽¹⁾(t)=IFFT[X(ω)]  (Expression 12)

In this expression, IFFT [ ] represents the inverse Fourier transform processing. The number of elements of a time signal vector x⁽¹⁾(t) obtained by the inverse Fourier transform is equal to the number of transmitting antennas. The (1) which is the superscript of x⁽¹⁾(t) indicates the first repetition process. The n^(th) repetition process is indicated as (n). In the following explanation, the superscript in other variables indicates the number of repetition processes.

The time signal vector x⁽¹⁾(t) obtained by the inverse Fourier transform 108 is connected to the memory 109 and the subtraction 110. The memory 109 stores the input signal x⁽¹⁾(t) only during the first repetition process. The output of the memory 109 is connected to the subtraction 110 and the selection 111. The processing of the subtraction 110 is necessary at the time of the second and subsequent repetition processes and will be described later.

In the selection 111, the signal from the memory 110 is selected and outputted when the repetition process is the first time. In the subsequent repetition processes, the signal from the subtraction 110 is selected and outputted to the Polar coordinate transformation 112.

In the polar coordinate transformation 112, the amplitude and phase of the time signal vector x⁽¹⁾(t) are calculated (Expression 13).

r ⁽¹⁾(t)=|x ⁽¹⁾(t)|

θ⁽¹⁾(t)=angle(x ⁽¹⁾(t))  (Expression 13)

It is necessary to perform the polar coordinate transformation processing for all signals of the respective elements and the respective time t.

The output of polar coordinate transformation 112 is inputted to the threshold value comparison 113 and the peak signal generation 114.

In the threshold value comparison 113, the amplitude r⁽¹⁾(t) is compared with a predetermined threshold value THR. If there is an amplitude exceeding the threshold value THR, it is determined that a large signal peak is present in the time signal of that time and element.

In the peak signal generation 114, a signal x_(p) ⁽¹⁾(t) indicating the peak component is calculated. The amplitude signal r⁽¹⁾(t) and the phase signal θ⁽¹⁾(t) from the polar coordinate transformation 112 are inputted to the peak signal generation 114. A peak signal vector x_(p) ⁽⁰⁾(t) represented by Expression 14 is calculated on the basis of the peak detection signal from the threshold value comparison 113.

$\begin{matrix} \left\{ \begin{matrix} {{x_{p}^{(1)}(t)} = {\left\{ {{r^{(1)}(t)} - {THR}} \right\} \cdot e^{j\; {\theta^{(1)}{(t)}}}}} & {{{if}\mspace{14mu} {r^{(1)}(t)}} > {THR}} \\ {{x_{p}^{(1)}(t)} = 0} & {otherwise} \end{matrix} \right. & \left( {{Expression}\mspace{14mu} 14} \right) \end{matrix}$

The {R⁽¹⁾(t)-THR} in the brackets in the upper equation of Expression 14 denotes the amplitude component of the peak exceeding a threshold value, and the e^(jθ) ⁽¹⁾ ^((t)) denotes the phase component thereof.

FIG. 4 is a conceptual diagram of Expression 14. The diagram shows a real number or an imaginary number component of the time signal vector x⁽¹⁾(t) and provides a threshold value THR is provided for the signal. However, the actual process shown in Expression 14 is performed on a complex number.

The peak signal x_(p) ⁽¹⁾(t) is inputted to the Fourier transform 115. In the Fourier transform 115, as in the case of the inverse Fourier transform 108, a Fourier transform is performed with a point number which is an integer multiple of KFFT. The peak frequency vector x_(p) ⁽¹⁾(ω) obtained by transforming the peak signal x_(p) ⁽¹⁾(t) into the frequency domain is calculated.

In the right singular value complex conjugate transpose matrix 116, the peak frequency vector x_(p) ⁽¹⁾(ω) is multiplied by the right singular value complex conjugate transpose matrix V^(H)(ω). As shown in FIG. 5, the stream vector S(ω) is transformed into the transmission signal vector X(ω) by multiplying the stream vector S(ω) by the right singular value matrix V(ω). Conversely, in the right singular value complex conjugate transpose matrix 116, by multiplying the right singular value complex conjugate transpose matrix V^(H)(ω), the peak frequency vector X_(p) ⁽¹⁾(ω) is decomposed into a stream region. A signal obtained by transforming the peak frequency vector X_(p) ⁽¹⁾(ω) into a stream region is represented by Expression 15. S_(p) ⁽¹⁾(ω) is referred to as a peak stream vector.

S _(p) ⁽¹⁾(ω)=V ^(H)(ω)X _(p) ⁽¹⁾(ω)  (Expression 15)

As described above, the peak stream vector S_(p) ⁽¹⁾(ω) is a signal obtained by transforming the peak signal x_(p) ⁽¹⁾(t) of the transmission signal vector x⁽¹⁾(t) into a frequency domain signal for each stream.

The peak stream vector S_(p) ⁽¹⁾(ω) from the right singular value complex conjugate transpose matrix 116 is inputted to the out-of-band level suppression 117. The peak stream vector S_(p) ⁽¹⁾(ω) is a function of the frequency ω. Since the peak stream vector S_(p) ⁽¹⁾(ω) is a wide band, it has a component even outside the band. Therefore, in the NULL carrier region outside the band shown in FIG. 6, if the value of the peak stream vector S_(p) ⁽¹⁾(ω) is replaced by 0, leakage does not occur outside the band.

The output of the out-of-band level suppression 117 is inputted to the stream level control integration 118. The amplitude level of the band-limited peak stream vector S_(p) ⁽¹⁾(ω) is controlled in accordance with the modulation multi-value number from the stream multi-value number control 102. The peak of the OFDM signal can be reduced by subtracting the peak signal decomposed for each stream from the signal not subjected to peak reduction. At that time, in order to prevent the overall bit error rate from being greatly deteriorated, processing is performed while appropriately adjusting the level thereof for each stream.

The detailed configuration of the stream level control integration 118 is shown in FIG. 2.

The peak stream vector S_(p) ⁽¹⁾(ω) from the out-of-band level suppression 117 is inputted to an adder 1181. In the adder 1181, the peak stream vector S_(p) ⁽¹⁾(ω) is added to the output signal from the memory 1184. In the first repetition process, the output result of the memory 1184 is 0, and the value of the peak stream vector S_(p) ⁽¹⁾(ω) is directly inputted to the clip 1183. The clip vector Q(ω)=[q₁(ω) q₂(ω) . . . q_(L)(ω)]^(T) from the clip value calculation 1182 is inputted to the other input of the clip 1183. In the clip 1183, the magnitude of each element of the clip vector Q(ω) is compared with the signal from the adder 1181. The magnitude of the addition result is limited so as not to exceed Q(ω). The method of calculating the magnitude of the clip value Q(ω) will be described later. The integration result I_(p) ⁽¹⁾(ω), which is the output of the stream level control integration 118, is finally equivalent to the value subtracted from the signal not subjected to peak reduction. The subtraction of the integration result corresponds to the application of a noise to a signal not subjected to the original peak reduction. Therefore, it is necessary to limit the amount of noise with the magnitude of the clip value Q(ω) in accordance with the multi-value number. For example, for a stream to which a large multi-value number such as 256 QAM or the like is allocated, the allowable EVM is also small. Thus, the value of the clip value q_(k)(ω) is also small. Conversely, in the case of a stream to which a small multi-value number such as QPSK or the like is allocated, the value of the clip value q_(k)(ω) may be large.

The output of the clip 1183 is outputted as the integration result I_(p) ⁽¹⁾(ω) and is inputted to the memory 1184. In the memory 1184, the integration result I_(p) ⁽¹⁾(ω) is stored in a storage element.

Brief description will be made on the second and subsequent repetition processes of the stream level control integration 118. The peak stream vector S_(p) ^((n))(ω) from the out-of-band level suppression 117 obtained in the n^(th) repetition process is inputted. In the adder 1181, the peak stream vector S_(p) ^((n))(ω) is added to the integration result I_(p) ^((n-1))(ω) of the previous cycle stored in the memory 1184. The clip 1183 and the memory 1184 perform the same process as the first repetition process and perform integration calculation of the peak stream vector S_(p) ^((n))(ω). The above process of the stream level control integration 118 is shown in Expression 16.

I _(p) ^((n))(ω)=[I _(p) ^((n-1))(ω)+S _(p) ^((n))(ω),Q(ω)]

I _(p) ⁽⁰⁾(ω)=0  (Expression 16)

In Expression 16, when the magnitude of each element of a vector A is greater than the magnitude of the corresponding element of vector Q(ω), B=clip[A, Q] represents the process of clipping with the element of the corresponding vector Q(ω). With respect to the processing of clip [A, Q], two specific examples will be described below.

As a first example, when a norm |a_(k)| of a_(k) is larger than q_(k), a_(k) is multiplied by a coefficient of q_(k)(ω)/|a_(k)| as shown in Expression 17. The level of the norm of a_(k) is suppressed so that the norm of a_(k) does not exceed q_(k).

$\begin{matrix} {b_{k} = \left\{ \begin{matrix} {a_{k} \times \frac{q_{k}(\omega)}{a_{k}}} & {{{if}\mspace{14mu} {a_{k}}} > {q_{k}(\omega)}} \\ a_{k} & {otherwise} \end{matrix} \right.} & \left( {{Expression}\mspace{14mu} 17} \right) \end{matrix}$

In this regard, A=[a₁ a₂ . . . A_(L)]^(T), and B=[b₁ b₂ . . . B_(L)]^(T). The elements of A and B are complex numbers, and the elements of Q are real numbers. However, in the processing of the first example, calculation and division of norms are necessary. There is a disadvantage that the calculation scale is large.

Therefore, the second example simplifying the calculation processing is shown in Expression 18.

$\begin{matrix} {{\left\lbrack b_{k} \right\rbrack} = \left\{ {{\begin{matrix} q_{k} & {{{if}\mspace{14mu} {\left\lbrack a_{k} \right\rbrack}} > q_{k}} \\ {- q_{k}} & {{{if}\mspace{14mu} {\left\lbrack a_{k} \right\rbrack}} < {- q_{k}}} \\ {\left\lbrack a_{k} \right\rbrack} & {otherwise} \end{matrix}{\left\lbrack b_{k} \right\rbrack}} = \left\{ \begin{matrix} q_{k} & {{{if}\mspace{14mu} {\left\lbrack a_{k} \right\rbrack}} > q_{k}} \\ {- q_{k}} & {{{if}\mspace{14mu} {\left\lbrack a_{k} \right\rbrack}} < {- q_{k}}} \\ {\left\lbrack a_{k} \right\rbrack} & {otherwise} \end{matrix} \right.} \right.} & \left( {{Expression}\mspace{14mu} 18} \right) \end{matrix}$

In Expression 18,

[ ] and ℑ[ ] are functions indicating a real number part and an imaginary number part.

In the calculation according to Expression 18, strictly speaking, the norm of a_(k) may exceed q_(k). However, the norm of a_(k) is close to q_(k), and the error is not a big problem. By the processing explained in Expression 17 or Expression 18, the output result of the clip 1183 is suppressed in level according to the modulation multi-value number.

The integration signal I_(p) ⁽¹⁾(ω) from the stream level control integration 118 is inputted to the selection 105.

The above process is the first repetition process. The process of decomposing the peak component of the OFDM signal for each frequency domain and each stream and controlling the level according to the multi-value number of the stream has been described above.

In the second and subsequent repetition processes, the selection 105 operates to output a signal from the stream level control integration 118. Accordingly, in the second and subsequent repetition processes, the output of the selection 105 is I_(p) ⁽¹⁾(ω), and the second and subsequent repetition processes are performed with respect to the peak stream signal.

The processing of the power distribution matrix 106 and the right singular value matrix 107 is similar to the first repetition process, and the output of the right singular value matrix 107 is represented by Expression 19.

X _(p) ⁽²⁾(ω)=V(ω)P(ω)I _(p) ⁽¹⁾(ω)  (Expression 19)

The peak frequency vector X_(p) ⁽²⁾(ω) is inputted to the inverse Fourier transform 108 and is transformed from the signal in the frequency domain into the signal in the time domain while performing oversampling as in the first repetition process. This inverse Fourier transform processing is shown in Expression 20.

x _(p) ⁽²⁾(t)−IFFT[X _(p) ⁽²⁾(ω)]  (Expression 20)

The peak time signal vector X_(p) ⁽²⁾(t) obtained by the inverse Fourier transform 108 is inputted to the subtraction 110. In the subtraction 110, as shown in Expression 21, the peak component is reduced by subtracting the peak time signal vector x_(p) ⁽²⁾(t) from the original time signal x⁽¹⁾(t) obtained in the first repetition process stored in the memory 109 and not subjected to peak reduction.

x ⁽²⁾(t)=x ⁽¹⁾(t) . . . x _(p) ⁽²⁾(t)  (Expression 21)

Although the peak time signal vector x_(p) ⁽²⁾(t) indicates the peak component of the time signal vector x⁽¹⁾(t), since the nonlinear calculation processing is performed in which the peak time signal vector x_(p) ⁽²⁾(t) is replaced by 0 in the out-of-band level suppression 117, it is not necessarily required that the peak signal can be reproduced accurately.

Therefore, even if the peak component is subtracted by Expression 21, the peak may remain although the peak component of x⁽²⁾(t) which is the subtraction result is reduced. Furthermore, by performing the nonlinear calculation processing, a peak component which has not been generated in the first repetition process may be newly generated. Therefore, it is impossible to completely reduce the peak by the process shown in Expression 21.

The transmission time signal vector x⁽²⁾(t) from the subtraction 110 is inputted to the selection 111. In the selection 111, the signal from the subtraction 111 is selected and outputted in the second and subsequent repetition processes.

In the polar coordinate transformation 112, the amplitude signal r⁽²⁾(t) and the phase signal θ⁽²⁾(t) of x⁽²⁾(t) are calculated. The amplitude signal r⁽²⁾(t) is compared with a threshold value THR in the threshold value comparison 113. The peak time signal vector x_(p) ⁽²⁾(t) is outputted in the peak signal generation 114. The peak time signal vector x_(p) ⁽²⁾(t) shows a peak component that cannot be reduced by Expression 21. The peak time signal vector x_(p) ⁽²⁾(t) is transformed into a signal in the frequency domain by the Fourier transform 115, and x_(p) ⁽²⁾(ω) is calculated.

The peak frequency vector X_(p) ⁽²⁾(ω) is again multiplied by the right singular value complex conjugate transpose matrix V^(H)(ω) in the right singular value complex conjugate transpose matrix 116 to output a peak stream vector S_(p) ⁽²⁾(ω).

In the out-of-band level suppression 117, the out-of-band component is replaced by 0 again. In the stream level control integration 118, integration and cropping shown in Expression 22 are performed on the first integration result I_(p) ⁽¹⁾(ω), and the result is outputted to the selection 105.

I _(p) ⁽²⁾(ω)=clip[I _(p) ⁽¹⁾(ω)+S _(p) ⁽²⁾(ω),Q]  (Expression 22)

The above processing is the second repetitive calculation process.

The third and subsequent repetitive calculation processes also perform the same calculation as the second repetitive calculation process. The peak component of the time signal x^((n))(t) is gradually reduced as the number of repetitive calculation increases. The peak time signal vector x_(p) ^((n))(t) and the peak frequency vector X_(p) ^((n))(ω) come close to zero. On the other hand, as the number of repetitive calculation increases, the peak component is integrated in the stream level control integration 118. The integration result I_(p) ^((n))(ω) converges to a value for reproducing the peak component.

By the series of processes, the peak component of the time signal x^((n))(t) is reduced, the signal amplitude reaching the nonlinear region of the power amplifier is reduced, and the nonlinear distortion component of the power amplifier is reduced.

Finally, a method of calculating the clip vector Q(ω) in the stream level control integration 118 omitted in the above description will be described. The peak stream vector S_(p) ^((n))(ω) is a signal obtained by decomposing the peak component into a stream. The peak stream integration vector I_(p) ^((n))(ω) obtained by integrating the peak stream vector S_(p) ^((n))(ω) with the number of repetition times is finally transformed into the peak signal x_(p) ^((n))(t) in the time domain. Then, by subtracting this signal from the time signal x⁽¹⁾(t) before reducing the peak, the peak component is reduced. As described above, the subtraction is performed in the time domain. Since this calculation is a linear process, as shown in FIG. 7, the calculation is equivalent to the subtraction of the peak stream integration vector I_(p) ^((N))(ω) from the stream vector S(ω) in the frequency domain decomposed into the stream.

As described above, the subtraction of the peak stream integration vector I_(p) ^((n))(ω) acts like a noise for the stream vector S(ω). A multi-value number is determined according to the magnitude of the eigenvalue Σ for the stream vector S(ω). The required SNR differs for each stream. For example, the required SNR of a stream using 256 QAM is about 17 dB, and the required SNR of a stream using QPSK is about dB. Therefore, if the peak stream integration vector I_(p) ^((n))(ω) of a large level is subtracted from the stream vector S(ω) to which 256 QAM is allocated, the EVM is greatly deteriorated, which leads to an increase in the bit error rate. On the other hand, if the stream has a low multi-value number such as QPSK or the like, even when subtracting the peak stream integration vector Ip(n) (ω) of a certain high level, no significant deterioration in the bit error rate occurs. In this manner, the magnitude of the peak stream integration vector I_(p) ^((n))(ω) that can be applied varies depending on the multi-value number allocated to the stream vector S(ω). The level thereof is suppressed by the clip vector Q(ω). As described above, in the stream level control integration 118, it is necessary to appropriately set the clip vector Q(ω) according to the modulation multi-value number. In the following, a method of calculating this clip vector Q(ω) will be described.

First, the reception SNR of the k^(th) stream will be mentioned below. Assuming that the power of the signal component in the reception signal is P_(signal) (k, ω), Expression 23 is established.

$\begin{matrix} \begin{matrix} {{P_{signal}\left( {k,\omega} \right)} = {{\lambda_{k}(\omega)}{P_{k}(\omega)}{E\left\lbrack {s_{k}}^{2} \right\rbrack}}} \\ {= {{\lambda_{k}(\omega)}{P_{k}(\omega)}P_{s}}} \end{matrix} & \left( {{Expression}\mspace{14mu} 23} \right) \\ {{P_{{peak}\mspace{14mu} {noise}}\left( {k,\omega} \right)} = {\sigma^{2} + {{\lambda_{k}(\omega)}{P_{k}(\omega)}{{i_{p,k}^{(n)}(\omega)}}^{2}}}} & \left( {{Expression}\mspace{14mu} 24} \right) \end{matrix}$

σ² in the first term of Expression 23 denotes the thermal noise power generated in the reception unit. The second term denotes the power of the k^(th) stream in the peak stream integration vector I_(p) ^((n))(ω). As described above, the peak stream integration vector I_(p) ^((n))(ω) acts as a noise. Thus, the peak stream integration vector I_(p) ^((n))(ω) is applied as a noise component on the transmission side and is a factor causing EVM deterioration.

From Expression 23 and Expression 24, the SNR of the k^(th) stream is represented by Expression 25.

$\begin{matrix} {{{{SNR}\left( {k,\omega} \right)} = \frac{P_{signal}\left( {k,\omega} \right)}{P_{{peak}\mspace{14mu} {noise}}\left( {k,\omega} \right)}}{{{SNR}\left( {k,\omega} \right)} = \frac{{\lambda_{k}(\omega)}{P_{k}(\omega)}P_{s}}{\sigma^{2} + {{\lambda_{k}(\omega)}{P_{k}(\omega)}{{i_{p,k}^{(n)}(\omega)}}^{2}}}}} & \left( {{Expression}\mspace{14mu} 25} \right) \end{matrix}$

In Expression 25, the required SNR not reducing the peak signal when |i_(p,k) ^((n))(ω)|²=0 is defined as Γ₀ and is represented by Expression 26.

$\begin{matrix} {{{SNR}_{0}\left( {k,\omega} \right)} = {\frac{{\lambda_{k}(\omega)}{P_{k}(\omega)}P_{s}}{\sigma_{0}^{2}} = \Gamma_{0}}} & \left( {{Expression}\mspace{14mu} 26} \right) \end{matrix}$

Further, the amount of deterioration of the required SNR due to peak reduction is assumed to be ε. The stream level control integration 118 operates so that the amount of deterioration of the required SNR falls below ε. For example, if the required SNR of 256 QAM is 17 dB and if the amount of deterioration of the required SNR due to peak reduction is 0.1 dB, 10 log₁₀(Γ₀)=17 dB and 10 log₁₀(ε)=0.1. Therefore, the amount of deterioration of the required SNR is defined by Expression 27.

10 log₀ {SNR ⁰(k,ω)}−10 log₁₀ {SNR(k,ω)}<10 log₁₀[ε]  (Expression 27)

Expression 27 is expanded and rearranged into Expression 28.

$\begin{matrix} {{{{10\; \log_{10}\left\{ \frac{{\lambda_{k}(\omega)}{P_{k}(\omega)}P_{s}}{\sigma_{0}^{2}} \right\}} - {10\; \log_{10}\left\{ \frac{{\lambda_{k}(\omega)}{P_{k}(\omega)}P_{s}}{\sigma_{0}^{2} + {{\lambda_{k}(\omega)}{P_{k}(\omega)}{{i_{p,k}^{(n)}(\omega)}}^{2}}} \right\}}} < ɛ}\mspace{20mu} {{10\; \log_{10}\left\{ \frac{\sigma_{0}^{2} + {{\lambda_{k}(\omega)}{P_{k}(\omega)}{{i_{p,k}^{(n)}(\omega)}}^{2}}}{\sigma_{0}^{2}} \right\}} < ɛ}\mspace{20mu} {{10\; \log_{10}\left\{ {1 + {\frac{\Gamma_{0}}{P_{s}}{{i_{p,k}^{(n)}(\omega)}}^{2}}} \right\}} < ɛ}\mspace{20mu} {{1 + {\frac{\Gamma_{0}}{P_{s}}{{i_{p,k}^{(n)}(\omega)}}^{2}}} < 10^{ɛ/10}}\mspace{20mu} {{{i_{p,k}^{(n)}(\omega)}}^{2} < {\frac{P_{s}}{\Gamma_{0}}\left( {10^{ɛ/10} - 1} \right)}}} & \left( {{Expression}\mspace{14mu} 28} \right) \end{matrix}$

In Expression 28, i_(p,k) ^((n))(ω) denotes the kth element of the peak stream integration vector I_(p) ^((n))(ω).

According to Equation 28, the level of the applicable peak stream integration vector I_(p) ^((n))(ω) is determined by the required SNR Γ0 and the allowable required SNR deterioration amount ε. Therefore, in the clip value calculation 1182, the clip value vector Q(ω) obtained according to Expression 28 is obtained. If the level is suppressed by the magnitude of the clip value vector Q(ω), the bit error rate is not greatly deteriorated.

$\begin{matrix} {{Q(\omega)} = {\frac{P_{s}}{\Gamma_{0}(\omega)}\left( {10^{ɛ/10} - 1} \right)}} & \left( {{Expression}\mspace{14mu} 29} \right) \end{matrix}$

FIG. 8 shows the summarized results of (10^(ε/10)−1)/Γ₀(ω) of Expression 29 required for calculating the required SNR: Γ₀ and the clip value vector Q(ω) for each modulation scheme of each stream on the premise of error correction at a coding rate 1/2, and examples of the allowable EVM deterioration amount.

In addition, the maximum number L_(max) of streams that can be transmitted is limited by the smaller value of the number of transmitting antennas N_(T) and the number of receiving antennas N_(R). Depending on the state of the transmission path characteristics H(ω), the number of transmitted streams L may be smaller than the L_(max). For example, when the number of antennas is four, the number of transmitted streams L becomes 3 or less in most of the transmission path characteristics H(ω). This is because in the fourth stream having the smallest eigenvalue, the eigenvalue becomes too small, whereby in most cases, it is difficult to perform transmission using this stream. In such a case, the fourth stream is not transmitted. Therefore, q4(ω) may be a large value or may satisfy 4(ω)≈Pusher 300. It is therefore possible to expect a large peak reduction effect.

According to the first embodiment of the present invention described above, by concentrating the component that reduces the peak of OFDM on the stream having a small eigenvalue, it is possible to effectively reduce the peak power of OFDM without largely deteriorating the bit error rate of the entire eigenmode MIMO transmission.

This application claims the benefit of priority based on Japanese Patent Application No. 2015-140473 filed on Jul. 14, 2015, the disclosure of which is incorporated herein by reference in its entirety.

INDUSTRIAL APPLICABILITY

The present invention can also be applied to an application for reduction of the peak power by concentrating a peak reducing component on a stream having a small eigenvalue.

Description of Reference Numerals  101: singular value decomposition  102: stream multi-value number control  103: mapping  104: 0 padding  105: selection  106: power distribution matrix  107: right singular value matrix  108: inverse Fourier transform  109: memory  110: subtraction  111: selection  112: polar coordinate transformation  113: threshold value comparison  114: peak signal generation  115: Fourier transform  116: right singular value complex conjugate transpose matrix  117: out-of-band level suppression  118: stream level integration 1181: adder 1182: clip value calculation 1183: clip 1184: memory 

What is claimed is:
 1. A peak power reduction device, comprising: a unit for dividing digital information to be transmitted into a plurality of streams; a unit for selecting a modulation multi-value number of the plurality of streams and distribution of transmission power according to a transmission state; a unit for singular-value-decomposing transmission path characteristics of the streams and precoding the transmission path characteristics of the streams by a right singular value matrix; a unit for complex-mapping subcarriers of the streams according to the modulation multi-value number; a unit for converting a complex mapping signal into a time domain signal; a first unit for storing a conversion result as the time domain signal; a second unit for, when a peak value exceeding a predetermined threshold exists with respect to an amplitude of the time domain signal, calculating a peak time signal exceeding a predetermined threshold value from the peak value; a third unit for converting the peak time signal into a frequency domain signal and calculating a converted frequency domain peak signal; a fourth unit for calculating a peak stream signal by multiplying the frequency domain peak signal by a complex conjugate transpose matrix of a right singular value matrix; and a fifth unit for limiting an SN ratio of a peak stream signal level corresponding to the modulation multi-value number of the streams with respect to the peak stream signal obtained by the fourth unit to a predetermined deterioration amount or less.
 2. The device of claim 1, further comprising: a sixth unit for calculating a peak time signal whose level is restricted by converting the level-restricted peak stream signal to a time domain; and a seventh unit for reducing a peak component by subtracting the time domain peak signal obtained by the sixth unit from the time domain signal obtained by the first unit.
 3. The device of claim 1, wherein a modulation error ratio for each stream is controlled according to the modulation multi-value number allocated to the streams by repeatedly operating the second unit to the seventh unit to calculate a peak-reduced time domain signal multiple times.
 4. The device of claim 2, wherein a modulation error ratio for each stream is controlled according to the modulation multi-value number allocated to the streams by repeatedly operating the second unit to the seventh unit to calculate a peak-reduced time domain signal multiple times. 